YITP-00-27 
June 2000 



On Electric Fields in Low Temperature Superconductors 



Jan GovaertsQ 

C.N. Yang Institute for Theoretical Physics 
State University of New York at Stony Brook 
Stony Brook NY 11794-3840, USA 

Damien Bertrand and Geoffrey Stenuit^ 

Institut de Physique Nucleaire, Universite catholique de Louvain 
2, Chemin du Cyclotron, B-1348 Louvain-la-Neuve, Belgium 



Abstract 

The manifestly Lorentz covariant Landau- Ginzburg equations coupled to Maxwell's equa- 
tions are considered as a possible framework for the effective description of the interactions 
between low temperature superconductors and magnetic as well as electric fields. A spe- 
cific experimental set-up, involving a nanoscopic superconductor and only static applied 
fields whose geometry is crucial however, is described, which should allow to confirm 
or invalidate the covariant model through the determination of the temperature depen- 
dency of the critical magnetic-electric field phase diagram and the identification of some 
distinctive features it should display. 
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1. Introduction. It is a widely held belief that (low T c ) superconductors cannot sustain 
electric fields in static configurations. The argument Q is directly based on the first of the 
London equations, 

^ (AJ) = E , dx (AJ) = —B , (1) 

where A is a phenomenological parameter proper to the superconducting material, J is the 
supercurrent density, and E, B are the electric and magnetic fields, respectively. Indeed, given 
that relation as well as the property of infinite conductivity, any nonvanishing electric field 
within the sample must set into motion a dissipationless supercurrent, hence a displacement 
of charges which very rapidly leads to an exact screening of any such electric field, certainly 
for time independent configurations (the case of stationary configurations with externally 
sustained supercurrents in the presence of magnetic vortices in type II superconductors is a 
different matter, see for example Refs.[§, ^ @). 

However, such a situation raises a series of puzzles of varying degrees of concern. First, 
it is physically inconceivable that an applied electric field would discontinuously drop to zero 
from its external value when crossing the surface of a superconducting sample. There ought 
to exist some skin effect with a characteristic nonvanishing penetration depth however small. 
Nevertheless, none of the parameters appearing in the London equations provides for such an 
electric penetration length. Moreover, in the above picture for the screening of electric fields, 
mention is made of an electromagnetic charge density which also is not accounted for in the 
London equations, nor more generally in the Landau-Ginzburg (LG) equations (see below). 

Another concern cit ci more formal level is the fact that the coupling of the London 
and LG equations to the electromagnetic fields is not spacetime covariant. Indeed, under 
Lorentz boosts, the supercurrent density J ought to transform as the space components of 
a 4-vector whose time component would then play the role of the aforementioned missing 
charge density, while the electric E and magnetic B fields transform as components of the 
two index antisymmetric field strength tensor F^ v . To illustrate this point more vividly 
perhaps, consider a flat infinite superconducting slab submitted to an external homogeneous 
magnetic field lying parallel to its surface. In such a case, the magnetic field will only partially 
penetrate the sample, with a characteristic penetration length related to the parameter A 
above Imagine now performing a Lorentz boost in a direction both parallel to the surface of 
the slab and perpendicular to the applied magnetic field. According to the Lorentz covariance 
of Maxwell's equations, in the boosted frame there appears now an electric field perpendicular 
to the surface of the sample, also within the volume of the superconductor where the magnetic 
held in the initial frame is nonvanishing. Thus Lorentz covariance requires the possibility of 
electric fields on the same footing as magnetic ones within superconductors, with an electric 
penetration length equal to the familiar magnetic one. Nevertheless, the existence of such 
electric fields within superconductors is incompatible with the London equations. 

One may take issue with the above covariance argument, since the superconducting 
sample itself defines a preferred frame, thereby insisting that the physics should be described 
only with respect to that specific rest-frame. Even though that frame is obviously distin- 
guished, the coupling of the London and LG equations to the electromagnetic fields should 
be consistent with the covariance properties of Maxwell's equations, even within that frame. 
In addition, one would also like to have available a manifestly covariant framework in which 
to study the interactions of electromagnetic fields and moving superconducting samples, an 
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issue which the usual London and LG equations are unable to address, and which is certainly 
accounted for to anyone's satisfaction in the description of electromagnetic interactions with 
ordinary conductors. It remains true nevertheless that some physical characterizations of 
superconductors can be defined with respect only to the rest-frame, such as for example 
the frame dependent notion of the free energy whose value determines the occurrence of the 
superconducting-normal phase transition only when evaluated in the rest-frame. 

Note that any spacetime covariant extension of the usual London and LG equations 
entails a time dependent LG (TDLG) equation in which space and time variations are on 
the same footing. Namely, a covariant TDLG equation is necessarily of second-order in time 
derivatives as it is in space derivatives, in sharp contrast with the usual non covariant first- 
order TDLG equations encountered in the literature jj]]. In particular, in a covariant setting, 
the time scale associated to time dependent fluctuations is then naturally set by the time it 
takes light to travel the distance of some mean value of the penetration and coherence lengths. 
In contradistinction for first-order TDLG equations, this time scale is specified in terms of 
an additional parameter, the relaxation constant |]]]. As a matter of fact, the latter quantity 
involves Boltzmann's constant, showing that the usual TDLG equations apply rather to the 
time dependency encured through thermodynamic fluctuations in the superconducting order 
parameter. As such, these TDLG equations do not provide a framework in which to study the 
intrinsically genuine time dependent dynamics of superconductors coupled to time varying 
electromagnetic fields, in the absence of thermodynamic fluctuations. A covariant extension 
of the London and LG equations would also provide such a dynamic framework, of relevance 
for instance to the dynamics of ensembles of magnetic vortices interacting among one another 
and with their electromagnetic environment. 

A manifestly covariant extension of the LG equation which immediately comes to one's 
mind is of course the so-called U(l) Higgs model of particle physics, whose construction itself 
was motivated by the BCS and LG theories in the late 50's. Indeed, as an effective theory for 
superconductivity, this model coincides with the original LG formulation for stationnary con- 
figurations, and readily provides || a description of all the remarkable quantum phenomena 
of superconductivity. It may thus appear somewhat surprising that the covariant formulation 
has not been used to further explore superconducting phenomena possibly lying beyond the 
boundaries of the usual London and LG equations. The purpose of this Letter is to suggest 
examples of investigations along such lines. 

After some considerations on the covariant London and LG equations presented in the 
next section, section 3 identifies a specific set-up which should enable to establish experi- 
mentally whether the covariant, rather than the usual noncovariant approach is relevant to 
superconducting phenomena in the presence of electric fields. What appears to us to be quite 
a remarkable circumstance is that this experimental confirmation of the covariant model 
should prove to be possible already using only a static, time independent configuration of 
external fields, through the observation of the superconducting-normal phase transition and 
the determination of the phase diagram in the (B,E) plane for a specific geometry of the 
applied fields. 

2. The covariant LG equations. As mentioned already, the considered model, of appli- 
cation only to low temperature superconductors, is that of a U(l) gauge invariant coupling 
to the electromagnetic interactions of a complex scalar field if) of charge q = — 2e < (the 
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Cooper pair charge) and with self-interactions determined through the usual LG potential 
(\ijj\ 2 — l) 2 properly normalized. The order parameter ip is normalized to the square-root of 
the Cooper pair density in the bulk in the absence of any electromagnetic field (see Ref.Q 
for some further details) . Rather than listing all the relevant equations in terms of the phys- 
ical quantities, let us already use the following choice of units. Space and time coordinates, 
namely x and the combination x° = ct with c being the speed of light in vacuum, are mea- 
sured in units of the penetration length A(T). Similarly, magnetic B and electric E/c fields 
are measured in units of <&o/(27rA 2 (T)), where $o = 27rft/|g| is the usual quantum of flux. 
Note that these units are temperature dependent, since the penetration length A(T) is, with a 
dependency we shall model]]} through A(T) = A(0) (l - (T/T c ) 4 ) 1/2 , T c b eing of course the 
critical temperature. Finally, the order parameter tp is parametrized according to ip = fe ld , 
with / real and f 2 thus measuring the relative Cooper pair density. In terms of these units, 
space and time coordinates are denoted u and r, and the magnetic and electric fields b and 
e, respectively. Finally, let us also introduce the quantities, 

A 3 (T) g\ 3 (T) 

h~F~ _/Wem ' J ~h~J 2 



. n q a {i ) -? q a~{± ) -> 

J = ^ , 2 PO CPem , 3 = T f2 J em , (2) 



where fio is the usual vacuum magnetic permitivity, and (cp em , J em ) are the superconducting 
electromagnetic charge and current densities (constructed in terms of ip), indeed defining a 
4- vector under Lorentz transformations. Note that these relations show that (/ 2 j°, f 2 j) is in 
fact proportional to this electromagnetic 4-supercurrent, which must be locally conserved. 

The latter remark is also confirmed by the inhomogeneous Maxwell equations (all space 
derivatives are of course with respect to u), 

d-e = -f 2 j° , dxb-d T e = -f 2 j , (3) 



which indeed, as is usual, imply the local conservation property d T (f 2 j°) + d ■ yf 2 jj = 0. 
The remaining electromagnetic equations of motion are given by, 

d T j + 8j° = -e , dxj = b, (4) 

which are recognized as the appropriate covariant extension of the London equations in (|l]). 
Note that only the first London equation is modified by the inclusion of a contribution of the 
supercharge density, as was indeed required, while the second London equation, essential to 
the Meissner effect, retains its original form. The homogeneous Maxwell equations, 

dxe + d T b = , d-b = 0, (5) 

follow from the covariant London equations (0) (as they do also from the noncovariant ones 
in ©). 

The equations of motion for the order parameter ip are given, on the one hand, by the 
covariant LG equation 



d 2 -d 2 \f = \f 2 -r\f- K \i-f 2 )f, ( 6 ) 

where the LG parameter k = A(T)/£(T) — £(T) being the coherence length — is essentially 
temperature independent|l|, and on the other hand, by the following conditions for the quan- 
tum phase 9, 

d T e = -j° + ip , dd = j-a. (7) 



3 



In these latter relations, <p and a are the scalar and vector gauge potentials defined such that 



e = —dip — d T a , b = dxa. (8) 

Finally, these equations are subject to boundary conditions requiring vanishing values for 
those components of the vectors df and j which are perpendicular to the surfaces of the 
superconducting sample in contact with an insulating material. 

The advantage of using this representation of the system is that the number of gauge 
dependent variables is kept to a minimum^]. Only the quantities 9, (p and a are defined up 
to the following gauge transformations, 

9' = 9 + X , tp' = ip + d TX , S' = a-d X , (9) 

x(u, t) being an arbitrary function, while the absolute sign of / may also be subject to gauge 
transformations [jTj . 

In fact, this decoupling of gauge variant and gauge invariant quantities may be ren- 
dered complete when substituting the covariant London equations (HI) into the inhomogeneous 
Maxwell equations In addition to the covariant LG equation (y), one then finds for the 

4-super current, 



d 2 - d 2 \ f = f 2 f - d T [d T f + d-j\ , [d 2 - d 2 \ j = f 2 j + d [d T f + d-j\ . (10) 

Any solution to this set of three coupled equations for j° , j and / then leads to specific values 
for b and e through the covariant London equations (|j), and in turn, once the gauge potentials 
(p and a related to these fields determined up to the gauge transformations parametrized by 
X, the corresponding solution for the quantum phase 9 is also finally obtained from (0)0. 

Although nonlinear, the equations @ and ([h]) also establish that these covariant LG 
equation admit progressive wave solutions with covariant dispersion relations in a linear 
regime. As a matter of fact, the phase and group velocities for variations in the 4-super cur rent 
and the order parameter / are different, unless the LG parameter k takes the same critical 
value k c = l/y/2 as the one which is so crucial to the understanding of the stability and 
interaction properties of magnetic vortices. Indeed, for a fluctuation of wave number k (in 
units of 1/| u\) around the vacuum solution (J° = 0, j = 0,f = 1), the group velocities are, 
respectively, 

Vj ~ 7PTT ' Vf ~ Vk 2 + 2k 2 ' 

thus showing that when k > k c (resp. k < k c ), waves in the 4-supercurrent will overtake (resp. 
be overtaken by) those in the order parameter /. In other words, within the superconductor, 
fluctuations in the electromagnetic fields will propagate more rapidly (resp. slowly) than 
those in the order parameter, in accordance with the relative magnetic and superconducting 
rigidities that the parameters A(T) and £(T) characterize. 

Clearly, such properties are totally different in the case of the usual noncovariant first- 
order TDLG equation, in which time scales are then normalized with respect to the relaxation 
parameter, which itself is temperature dependent. The ensuing dispersion relations are then 
linear in frequency, implying that the group velocities of fluctuations in the supercurrent j 
and in the order parameter / are then also identical, independently of the value for k. Such 
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differences between the covariant and noncovariant frameworks must lead to distinct physical 
properties in the case of time dependent configurations in ultra-high frequency regimes, v ~ 
c/A(T), c/£(T), an issue which, however, is beyond the scope of this work. 

To conclude this general discussion, let us also give the expression for the free energy 
E of the system in the covariant form, 

yy T e - l ^ + p- mi + 

+ l i'afaff + (3ff + f (f 2 + ?) + \^ (1 - f) 2 - , (i2) 

where the normalization factor related to our choice of units is displayed together with E in 
the l.h.s., e cx t and 6 ex t are externally applied electric and magnetic fields, respectively, and 
Q stands for the volume of the superconducting sample. 

The same expression is also of application to the noncovariant model, in which case 
one has j° = and e = within the superconductor, and the quadratic term in d T f is to 
be replaced by a linear term while the time coordinate is then also measured in units of the 
relaxation parameter for the TDLG equation. The term in [e — e cx t] 2 measures the energy 
required to expulse the electric field from the superconductor. In the noncovariant case, 
and in accordance with the first London equation, we shall thus assume that the associated 
penetration depth is essentially vanishing for all practical purposes. For physics reasons, 
such an approximation cannot be very reliable when it comes to nanoscopic superconductors, 
but we shall use it as a working hypothesis anyway. Note that the free energy E is defined 
here in such a way that it vanishes at the superconducting-normal phase transition. And 
as a last remark, clearly, in the case of stationary configurations and in the absence of any 
electric fields, the equations of both approaches coincide with the usual noncovariant time 
independent LG equations. 

3. Characterizing the phase transition. In order to identify a specific geometry of 
applied fields which could help discriminate experimentally between the two approaches al- 
ready in a static configuration, consider again the situation of the flat infinite slab of the 
Introduction, this time subjected not only to the homogeneous magnetic field parallel to its 
surface, but also to an homogeneous electric field applied perpendicularly to its surface (an 
electric field parallel to the slab does not induce a supercharge distribution j°, hence neither 
a feature distinctive from the noncovariant model). The slab is taken to be of thickness 2a, 
while the external electric field e ex t is aligned along the x axis, and the external magnetic 
field 6 cx t along the y axis, with components e ex t and b ex t, respectively (the origin of this co- 
ordinate system is of course positioned in the center of the slab). For this specific geometry, 
the expulsion of the magnetic field is achieved through an induced supercurrent j circulating 
along the z axis, also parallel to the slab, while that of the electric field is achieved through 
the appearance of a nonvanishing supercharge density j°, an occurrence which simply cannot 
arise in the noncovariant approach. Both these effects imply a deviation from its canonical 
value of unity for the order parameter /. In view of the symmetries of the problem, both 
j°(u) and j z (u) are odd functions of the normalized u = x/X(T) coordinate along the x axis, 
while f(u) is even. 
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Given the equations and the different conditions imposed on these quantities at the 
boundaries u = ±u a = ±a/A(T), it proves possible as well as useful to express both j°(u) 
and j z (u) in terms of a single function j{u) 



3° in) 



-e ex tJW 



(13) 



so that one has for the electric and magnetic fields within the sample, 



e(u) = e ext ^-j(u) , b(u) 
du 



du y ' ' 



(14) 



thus showing once again that Lorentz covariance implies that the penetration lengths for both 
types of fields are identical. The set of equations to be considered then reduces to, 



-/*<«»<«) , £/<«)-(< 

subject to the boundary conditions 



2 

ext 



j 2 (u)f(u)- K 2 (l-f(u))f(u) , (15) 



d 



du' 



. 



(16) 



Note already the subtle interplay between the magnetic and electric field contributions to the 
LG equation for f(u), which leads to values larger than unity for f(u) in the electric regime 
e cxt > fr<Lt> while, as is usual, f(u) remains less than unity in the magnetic regime b 2 xt > e 2 xt . 
The existence of these two regimes is a direct and distinctive consequence of manifest Lorentz 
covariance; only the magnetic one arises in the noncovariant approach (see below). Note also 
that in view of these equations, the solutions for j(u) and f(u) are necessarily functions of 

only, indeed justifying this notion of electric or magnetic 



ext 



-'ext. 



the specific combination (6, 
regimes. 

Up to the normalisation factor displayed in (12) as well as the infinite surface of the 
slab, the free energy £ of configurations obeying these equations is simply given by 



£ = 2u n 



1 j(u a ) 



Un 



f^ext + e ext) 



— / du 

U a JO 



2 

ext 



cxl; I .i~f ~ + l^^f 



(17) 



Consequently, the curve in the (6, e) phase diagram which characterizes the superconducting- 
normal phase transition obeys the following equation in the covariant approach, 



b 2 + e 2 



1 



i 

U a JO 



du 



(> 2 



(18) 



The noteworthy property of this relation is that the l.h.s. involves only the combination 
(b 2 + e 2 ), while the r.h.s. is only a function of the combination (b 2 — e 2 ) of the external fields 
(since the solutions j(u) and f(u) also share that property). 

Before addressing the specific consequences of this equation for the (b, e) phase diagram, 
let us consider the corresponding expressions in the noncovariant approach. In that case, one 
has of course j°(u) = (which also implies e(u) = within the superconductor, as follows 
from the first London equation) as well as j z (u) = —b cxt j(u). The equations then remain as 
given in (|l5|) , including the boundary conditions, with the only but important difference that 
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the factor (b 2 xt — e 2 xt ) appearing in the LG equation is of course replaced by 6 2 xt only. Hence, 
the noncovariant LG equations only admit the magnetic regime of solutions. Note that in 
this case, the solutions for j(u) and f(u) are then also functions of the 6 2 xt external field only. 
Consideration of the expression for the free energy then leads to the following condition of 
criticality in the (6, e) phase diagram in the noncovariant case, 



6 2 + 



1 



1 



i 



/ </?/ 



b 2 J 2 f + \, 2 f 



(19) 



In spite of the apparent similarity with (18), recall however that the r.h.s. of this expression is 
a function of b 2 only, while the l.h.s. is no longer the simple combination b 2 + e 2 characteristic 
of a circle since the coefficient multiplying the term in e 2 is also a function of b 2 . 

Note that the conditions (]l8| ) and ([H]) coincide in the limit that no electric field is 
applied, e = 0, as they should of course. Moreover in the absence of any magnetic field, 
6 = 0, ([ll]) implies the existence of a nonvanishing critical electric field, eo = k/V2, or 
in physical units E (T)/c = (A(0)/A(T)) 2 B™(0), B™(0) = $ /(2\/2vrA(0)^(0)) being the 
usual thermodynamic critical magnetic field in the bulk at zero temperature. Clearly, the 
existence of such a critical electric field even in the noncovariant approach is consequence 
of our definition for the free energy in ( |l2| ) which accounts for the expulsed electric energy 
density through the term in [e — e ex t] 2 . In particular, this critical electric field Eq(T) vanishes 
at the critical temperature T c , as does the critical magnetic field Bq(T) in the absence of any 
electric field, e = 0. 



4. The (B, E) phase diagram. A complete unravelling of the consequences of the criticality 
conditions (|l8| ) and ( |l9|) requires of course a numerical approach. Nevertheless, an analysis 
in some limiting situations already suffices to gain insight into the differences implied by the 
two models. An obvious such situation is obtained in the macroscopic limit, namely when the 
slab half-thickness a is much larger than both the penetration and coherence lengths. For all 
practical purposes, the function j(u) then essentially vanishes whereas the order parameter 
retains its canonical value of unity within most of the volume of the sample, except for a 
small region close to the surface. Hence in the above expressions of criticality, in the limit 
that a — > co, only the contribution in k 2 / 4 /2 tends to dominate, leading in both cases to the 
condition, 

b 2 + e 2 ~^K 2 , a»A(T),£(T). (20) 

Since this will prove to be useful, let us normalize the measurement of these fields to the 
value 6o of the critical magnetic field in the absence of any electric field (in the macroscopic 
limit, we thus have 6o = k/^/2). In terms of the physical quantities, one then obtains the 
following approximation to the criticality condition in the (B, E) phase diagram 

(|) 2 + (lf) 2 = 1 • a » A <^«^ < 21 > 

Hence in the macroscopic limit, the two models are not distinguished in their (B,E) phase 
diagrams. In particular, both their critical magnetic, Bq, and electric, Eq, fields (in the 
absence each time of the other field) reach a vanishing value at the critical temperature T c . 
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Consider now the nanoscopic limit, namely when a <C A(T),£(T). In practice, this 
situation may be encountered indeed for nanoscopic samples close to the critical temperature 
T c . In such a case, one may develop series expansion solutions in u for the functions j(u) and 
f(u) in order to evaluate the criticality conditions ( |l8|) and (|l9|). However, since contributions 
of order b 2 and e 2 appear on both sides of these equations, in order to be of sufficient accuracy, 
the expansion in u must at the same time include at least the first order corrections in b 2 
and e 2 in the r.h.s. of ( |l8| ) and ( |i~9|) as well. For this reason, it is more relevant to consider 
a weak field expansion for the solutions independently of whether u a is small or not, to be 
used to compute to first order in b 2 and e 2 the r.h.s. of the criticality conditions above, and 
then eventually take the nanoscopic limit. Note that given the result (^), critical fields are 
at least of the order of n/y/2, so that such a weak field expansion should be warranted for 
small values of the LG parameter n, namely for type I superconductors. 

After some work, one then finds that the criticality conditions (|l^) and (|i~9|), evaluated 
to first order in b 2 and e 2 , imply the following constraint on the physical fields in the (B, E) 
phase diagram, 



.Bo J \ B 

where as before Bq stands for the critical magnetic field value in the absence of any electric 
field, E = 0, which is in general a function of temperature and of a of course, while C is a 
factor given by the following expressions, 

covariant model : C = f > 

(23) 

noncovariant model : C = ( Ua _^ nhlla ) (j^) , 

where 

u„ 1 [ i— tanh 2 u n A 9 

13 = Tr? r-T \2 ( 2 ^2 W2 a - (3k 4 - 10« 2 + 16) tanh Ua + 

+ (5k 4 - 22k 2 + 16) tanh 3 u a + (k 2 - 2) (3k 2 - 4) — — — } . (24) 

cosh u a J 

These expressions are valid in the weak field approximation to first order whatever 
the value for a. Taking now the nanoscopic limit as well, one finds (1 — tanh(u a )/u a ) 1 = 
3/u 2 [l + 0(u 2 )] and (3 = 1/2 [l + 0(u 2 )], leading finally to the following criticality condi- 
tions in the (B,E) phase diagram in the weak field limit, 

covariant model : (^) 2 + 3 (^f - 1 , a < A(T),£(T) , 

noncovariant model : )' +6 (^ff —f—, i^gf ~ 1 , a « A(T)^(T) . 

\T C ) 

(25) 

Since the critical magnetic field Bo does vanish towards the critical temperature T = T c , so 
do all the critical fields B and E which are defined by either of these relations, and thus in 
particular also the critical electric field Eq in the absence of any magnetic field, B = 0, as 
was already remarked previously in the noncovariant case. However, by having chosen to 
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normalize the measurements of these fields to Bq, a very distinctive feature appears for the 
covariant model when compared to the noncovariant one. Indeed, the ratio E/(cBq) always 
retains a finite and nonvanishing value, whatever the critical values for B and E within 
the intervals [0, Bq] and [0, Eq], even in the limit of the critical temperature T c , whereas in 
the noncovariant model, that same ratio E/(cBq) must vanish like \A — {T /T C ) A (given our 
chosen model for A(T)). In particular, in the weak field approximation and including the 
result (^0|) valid for macroscopic samples, one thus derives in the covariant case the following 
bounds for the critical electric field Eq, 



with Eq/(cBq) moving towards lower values within that interval when the critical temperature 
T c is approached (recall that (3 is also temperature dependent through u a ). In the nanoscopic 
limit a <C A(T),£(T), these same bounds reduce to 

In contradistinction in the noncovariant case, the lower bound on Eq/(cBq) always vanishes, 
since one then finds, 

i/fl--tanhn a ) (1-/3) < ^ < 1 , (28) 

V V U a J Bq 

reducing in the nanoscopic limit a <C A(T), £(T) to 

V6U(o)/r vt c /<^ £ < 1 - < 29 > 

As a matter of fact, this type of consideration may be refined further still in the covari- 
ant case. Indeed, an obvious solution to the covariant LG equations is j(u) = sinhu/ coshu a , 
f(u) = 1 in the case that e ex t = ^ext, a fact which, as was remarked previously, is a dis- 
tinctive feature of the covariant approach, since this solution defines precisely the boundary 
between the magnetic and electric regimes of superconductivity, and as such its existence is 
a direct consequence of Lorentz covariance. Hence, the critical condition ( Jl8| ) for the corre- 
sponding fields bi and e\ simplifies in this specific instance to the exact result, valid under 
all circumstances, 

b\ + e\ = 2b\ = 2e\ = -k 2 ^ . (30) 

2 1 tanh u a 

When the weak field approximation is also warranted for the evaluation of the critical mag- 
netic field Bq, this result combines with those above to lead to the following bounds, 

1 , Bi E x lc 1 , . 




and in the nanoscopic limit, 

1 

2 < Bq ~ B " ^2 



l < |i = ^/c < 1 > a « A(TU(T)) (32) 
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whereas in the noncovariant case, one finds similarly in the nanoscopic limit 



<^ = ^<4 , a«A(T),C(T). (33) 



1 + 6 V • I 



/A(0)\ 2 1 B Bq y/2 



Hence here again for those specific configurations such that £> ext = E ex ^/c, the lower bound 
on Bi/Bq reaches a vanishing value at the critical temperature in the noncovariant case, 
whereas that lower bound remains finite and is only mildly temperature dependent in the 
covariant case. 

The existence of such finite bounds on the values for E\/(cB\) as a function of tem- 
perature in the covariant case, translates into the following nice characterization in terms of 
the (B/Bq,E/(cBq) phase diagram. Indeed, the limits (|32|) (which are more refined in pi])) 
imply that the phase boundary curve in that diagram must always cross the diagonal line 
B = E/c within the interval of B/Bq or E/(cBq) values defined by these bounds in the co- 
variant model, whatever the value for the temperature (see Fig.l). Such a property is simply 



not met in the noncovariant model (see Fig. 2). Similarly, the lower bounds (26) or (|27| ) on 
Eq/(cBq) imply that, when approaching the critical temperature, the same phase boundary 
curve at B/Bq = cannot move below a specific finite value in the covariant model, while it 
must necessarily do so in the noncovariant one. 

As a conclusion thus, which should remain valid beyond the specific limits conside- 
red here, it appears that by choosing to normalize the measurement of critical electric and 
magnetic fields to the critical magnetic field in the absence of any electric field, for a given 
nanoscopic sample with this specific geometry of applied fields and by approaching the critical 
temperature, the (B,E) phase diagram provides the necessary distinctive features which 
should enable to discriminate experimentally between the covariant and noncovariant mo- 
dels, and in any case confirm or invalidate the description offered by the Lorentz covariant 
LG equations. Indeed, as was remarked previously, the ordinary noncovariant framework is 
not physically realistic when it comes to nanoscopic samples in the presence of electric fields, 
since it ignores the partial penetration, albeit small, of the electric field into the sample's 
surface. The present analysis has concentrated on the weak field approximation, essentially 
in the nanoscopic limit. Similar distinctive differences between the two models should also 
exist for larger values of n, in ways still to be investigated requiring then a detailed numerical 
study which is not pursued in this Letter. 



5. Numerical solutions. Here, we present the results of the numerical resolution of the 
LG equations and of the criticality conditions (|l^) and (|l9| ) for only one situation, which is 
close enough both to the discussion of the previous section and to an experimentally realistic 
situation. Namely, we take the following parameter values 

m =5 • ,i=0 - 02 - <34) 

Indeed, this value for k is typical for aluminium (Al), while tabulated values of A(0) for 
Al — A(0) = 16 — 50 nm with T c = 1.18 K — would imply that the slab is then a few hundred 
nanometers thick, within reach of present lithographic techniques for Al on a Si02 substrate. 
Moreover, the critical magnetic field 2?£°(0) for Al is also on the order of 100 Gauss, so 
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that the required electric field values for a measurement of the (B, E) phase diagram would 
reach into 3 MV/m, namely 3 V//nm, certainly also a reasonable range of values for such 
a nanoscopic device. Of course, compared to the infinite slab model, such a device will be 
subjected to finite size corrections. Presumably, such corrections would imply that the role 
played by A(T) and k in our analysis would be replaced by some effective quantities whose 
values would not differ to a great extent from those of Al in the bulk. Such corrections may 
be assessed only once a specific device is designed. 

In Fig.l (resp. Fig. 2), we present the (B/Bq, E/(cBq)) phase diagram for the covariant 
(resp. noncovariant) model, given the values in (|34|), for a series of temperatures in the 
range from T = to T = T c . The general behaviour of the phase diagram as a function 
of temperature is indeed the one described in the previous section. In particular in the 
covariant model, and as a function of temperature, the critical electric field values Eq/{cBq) 
and E\/{cBq) obey the different finite lower (and upper bounds) derived from the analytical 
discussion, including those given in (26) and (|3l|) when considering the associated values for 
(3. In contradistinction, in the noncovariant case, the ratio E/(cBq) reaches a vanishing value 
when approaching the critical temperature, while in this case as well it way be checked that 
the different lower bounds (pl|) , (|29|) and (|33|) are indeed also obeyed. 

Such results, as well as the other considerations of this Letter show that it should be 
possible to experimentally discriminate between the ordinary noncovariant LG equations and 
the covariant ones advocated here, by determining the critical (B, E) phase diagram of a 
nanoscopic superconducting sample for temperatures approaching its critical temperature. 
The geometry of the applied fields is crucial for this purpose, with the external magnetic 
field parallel to the sample's surface and the external electric field perpendicular to it. By 
normalizing the measurement of fields to that of the critical magnetic field in the absence 
of any electric field, distinctive differences between the two approaches are best brought to 
the fore, and should enable to confirm or invalidate the covariant approach. Moreover, if the 
experiment should also allow for an absolute calibration of the applied fields, the comparison 
between the two models may be refined still further by considering the temperature depen- 
dency of the critical value for applied magnetic B and electric E/c fields of equal magnitude, 
this temperature dependency being constrained to lie within a specific interval whose exis- 
tence is a direct consequence of the manifest Lorentz covariance of the covariant model. We 
hope to be able to report on such measurements in the future, but lithographic problems 
have hindered any progress until now. 
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Figure Captions 



Figure 1: The phase diagram (B/Bq, E/(cBq)) for the covariant LG equations with the values 
(pip. Shown from top to bottom are the curves associated to the following increasing tem- 
perature values, T/T c = 0,0.8766,0.9659,0.9935,0.9996. The diagonal line determines those 
configurations such that B ex t = E ext /c, the two vertical dot-dashed lines at B/Bq = 1/2 and 
B/Bq = l/s/2 correspond to the lower and upper bounds ( |32|) obeyed by the critical electric 
Ei/c and magnetic B\ fields of equal strength in the nanoscopic limit of the weak field appro- 
ximation for the covariant LG equations, while the horizontal dashed line at E/(cBq) = l/y/3 
corresponds to the lower bound (^7j) on the critical electric field Eq/(cBq) in the same ap- 
proximation. The existence of these finite bounds is the distinctive prediction of the covariant 
model and a direct consequence of its manifest Lorentz covariance. 

Figure 2: The same as in Fig.l for the noncovariant model. In this case, the horizontal and 
two vertical lines are displayed only for the purpose of comparison with the covariant model. 



13 



Figure 1 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

B/BO 



Figure 2 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

B/BO 



